Spontaneous squeezing of a vortex in an optical lattice 
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We study the equilibrium states of a vortex in a Bose-Einstein condensate in a one-dimensional 
optical lattice. We find that quantum effects can be important and that it is even possible for the 
vortex to be strongly squeezed, which reflects itself in a different quantum mechanical uncertainty of 
the vortex position in two orthogonal directions. The latter is observable by measuring the atomic 
density after an expansion of the Bose-Einstein condensate in the lattice. 
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I. INTRODUCTION 

Vortices play a crucial role in explaining rotational and dissipative properties of supcrfluids. In the superfluid 
systems studied until now, such as liquid Helium 1], superconductors Q, and Bose-Einstein condensates 3, 4, 5, f| 
ItI IM lol lift fill ■ a vortex line behaves as a classical object. However, placing a Bose-Einstein condensate in an optical 
lattice |l2l Il3l IT3. ITsl ITij leads to an unprecedented control over the system parameters, which has already enabled 
experimental studies on quantum phase-transitions ^ISj 0, E|) superfluidity 12, 12j, and number squeezing .1 j |. 
Furthermore, a Bose-Einstein condensate in a one-dimensional optical lattice has a similar layered structure as the 
cuprate superconductors and is a promising way to achieve the quantum Hall regime in a Bose-Einstein condensed 
gas |l9j| . This suggests that a careful examination of the quantum properties of a vortex line in a one-dimensional 
optical lattice is warranted. 

A characteristic feature of all superfluids is their ability to support persistent currents. Quantized vortices play 
a crucial role in understanding the decay of these currents, and also the superfluid's response to rotation or to an 
external magnetic field in the case of superconductor. In pioneering work Fetter developed the quantum theory of 
vortices in liquid Helium 20] and so-called type-II superconductors l2l| . In these systems quantum fluctuations turned 
out to play such a small role that they are experimentally inaccessible. Here we show, however, that truly quantum 
mechanical behavior of the macroscopic vortex occurs in a Bose-Einstein condensate in a one-dimensional optical 
lattice. This is due to the reduced dimensionality as well as the reduced number of atoms in every well of the optical 
lattice. Remarkably, it turns out that the vortex can spontaneously become strongly squeezed, which is reflected in 
the quantum-mechanical probability distribution of the vortex position. Unlike coherent states, which are described 
by Gross-Pitaevskii theory, squeezed states are highly nonclassical. 

II. BOSE-HUBBARD HAMILTONIAN 

We consider a Bose-Einstein condensate in a one-dimensional optical lattice with lattice spacing A/2, where A is the 
wavelength of the laser beams creating the standing wave of the lattice. The optical lattice splits the Bose-Einstein 
condensate into a stack of Ns weakly-coupled pancake condensates, each containing N atoms. For concreteness we 
use in the following always ^^Rb atoms. The depth of the lattice Vl can be easily changed and controls the tunneling 
of an atom from one pancake condensate to the next, i.e., the strength of the interlayer Josephson coupling. The 
condensate also experiences a harmonic trapping potential in the radial direction with an oscillator frequency LOr- 
• ^ i The longitudinal trapping potential along the direction of the laser beams is assumed to be so weak that it can be 
neglected. While the lattice is taken to be deep enough to allow us to use a tight-binding approximation and to include 
^ , only the weak nearest-neighbor Josephson coupling, it is also taken to be shallow enough to support a superfluid state 
■ - - ' as opposed to the Mott-insulator state . 

Furthermore, we consider a vortex line that pierces at the position i/n) through each layer of the stack labelled 
by the index n. In each layer the density of the pancake condensate varies as a function of the radial distance. As a 
result the vortex experiences an effective potential that depends on its distance from the origin ■ Without rotation 
of the gas this potential is approximately an inverted parabola. Hence, the vortex is energetically unstable and will 
tend to spiral out of the system. However, the vortex is stabilized if we rotate the Bose-Einstein condensate with 
a rotation frequency Q that is larger than a critical frequency flc. In the rest of this paper we consider the regime 
where the vortex is energetically stable and Q > D,c. 
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In the tight-binding approximation the attraction between the nearest-neighbor parts of the vortex turns out to 
be harmonic with respect to their separation. We denote the typical strength of this attraction by Jy. Physically 
the attractive interaction is due to the energy cost for phase differences between the two layers. The theory of the 
resulting coupled harmonic oscillators can be quantized by introducing the bosonic annihilation operators for the 
eigenmodes of the vortex line 23^ 2M- These modes are the Kelvin modes (or kclvons) and correspond physically 
to a wiggling of the vortex line with a wavelength 27r/fc. The operators Vk in momentum space are related to the 
coordinate space kelvon operators and to the vortex positions through 

and 

yn = (2) 



2VN 

where R is the typical radial size of each pancake condensate. In this way we obtain the dispersion relation hLUpcik) = 
hujK{0) + Jy [1 — cos (fcA/2)] of the kelvons "2^. The energetic stability of the vortex is reflected in a positive value 
of the dispersion at zero momentum. Physically the frequency at zero momentum is the precession frequency of a 
slightly displaced straight vortex line around the center of the Bose-Einstein condensate. 

It is crucial for our purposes to realize that, due to the Euler dynamics of the vortex, the coordinates x„ and t/n 
are canonically conjugate variables and therefore obey the Hcisenberg uncertainty relation [a;„,y„] = iR^/2N. This 
means that the position of the vortex is always smeared out by quantum fluctuations. For a vortex in a lattice these 
fluctuations can become large due to the reduced particle numbers in every layer, as opposed to the total number of 
particles, and the radial spreading of the pancake condensate wave function as the lattice depth is increased. 

The kelvon dispersion relation was determined by taking into account only the lowest order expansion of the effective 
potential experienced by the vortex. However, expanding up to fourth order in the vortex displacements results in an 
interaction (Vb/2) v^vl^VnVn for the kelvons. In the regime where the vortex is stable Vq < and this interaction is 
attractive. The physics of the vortex line is thus described by a one-dimensional Bose-Hubbard model with a negative 
interaction strength. The corresponding Hamiltonian is given by 

H = [TiujKi^) + Jv] vivn - ^ XI Y ^ ^'i^i'^"'^"' (3) 

n {n,m) n 

where (n, m) denotes the nearest-neighbor layers. 

We calculate the parameters of this Bose-Hubbard model using a variational ansatz we used in our earlier work [23| , 
where the density of the condensate wave function in each layer was proportional to exp [— [x^ -f j/^) / ■ In this 
way we find that the strength of the nearest-neighbor coupling is given by 



Jv — -TT^ 



1 exp 



/ V 4 / V Ah 



and the interaction strength is given by 

2huJr (Ir/Rf r [0, If/R^] - ShcJr (Ir/Rf " 4?iO 



Vo 



where — \/ Stt'^Vl / mX^ is the oscillator frequency of the optical lattice, Ir — ^JhJrmTr^ and F [a, z\ is the incomplete 
gamma function |2!tI |. Furthermore, a straight and slightly displaced vortex precesses around the condensate center 
with the frequency 

ujKik = 0) = {cJrll/2R^) (1 - F [0, l^/R^] ) + Q. (4) 

The precession frequency of the vortex changes from negative to positive at the critical rotation frequency 

= {L0rfj2R^) (F [0, It/R"] - 1) . (5) 

When n > ilc the vortex is locally energetically stable. The condensate size R is given by the solution of the transcen- 
dental equation (Ir/R)^ 1 — 7 -I- 2Na^J7nIoZp2Tih — 4 In (Ir/R) = 1, where 7 is the Euler-Mascheroni constant and a 

is the three-dimensional scattering length. Without the kelvon interaction the ground state of the vortex corresponds 
to a vortex line along the symmetry axis of the Bose-Einstein condensate and with a quantum mechanical uncertainty 
such that (2;^) = (y^j). An attractive interaction will lead to squeezing with (x^) ^ (7/^) as we show next. 
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III. VORTEX SQUEEZING 

First we have to identify the correct order parameter for the vortex hne. Since H. > ilc there occurs no Bose-Einstein 
condensation of the kelvons and the appropriate order parameter for the vortex hne is not (vn), which would signal 
the above mentioned energetic instability and the tendency of the vortex to move away from the symmetry axis. 
However, the correct order parameter can be identified with A = Vo{vnVn)- To arrive at the associated mean-field 
theory, we quadratically expand the Hamiltonian around this order parameter and obtain apart from a constant 



H = Y, [Tiu;K{0) + Jv] vivn - ^ E + E 

n {n,m) n 



(6) 



We diagonalize this result by means of a Bogoliubov transformation. It is most convenient to work in the momentum 
space where we define new creation and annihilation operators with the help of of the Bogoliubov amplitudes Uk and 
Vk as bk = Ukvl — VkV-k and bl = u^Vk — vlv^_i^. Furthermore, to ensure bosonic commutation relations we must have 
|ufeP ^ l^'fcP = 1- The Hamiltonian in Eq. © is then diagonalized with the choise 



(£;(fc)+^K(fc)r-|AP 

and 



vk\'^ . ' (7) 



lu.l^ = , (8) 



where E{k) = y/ajKik)"^ — |Ap is the dispersion of the Bogoliubov quasiparticles. 
Requiring self-consistency of the approach leads to a gap equation 

Vo 2E{k) ' 

where Nk = \/{e^^^'^^ — 1) is the Bose distribution. We solve the gap equation for the order parameter, analogous 
to the BCS theory of superconductors jl^l . In Fig. ^ we present the behavior of the order parameter as a function 
of temperature and rotation frequency for typical parameter values. We see from Fig. ^ that there is a critical 
temperature below which the order parameter becomes nonzero and then increases monotonically with decreasing 
temperature. Furthermore, the absolute value of the order parameter at zero temperature is usually close to the 
precession frequency. Only very close to the critical rotation frequency Sic the deviation becomes large. 

The complex order parameter A — |A|e''^ has an interesting physical interpretation in terms of the quantum 
mechanical uncertainty of the vortex position. It turns out that in a coordinate system rotated by an angle 9 we have 



where 4> is the phase of the order parameter. This implies that a nonzero value of the order parameter is refiected 
in the squeezing of the vortex position distribution. The main axis of the uncertainty ellipse is at an angle = 0/2 
with the x-axis. In equilibrium the uncertainty ellipse of the vortex position distribution is independent of the layer 
index n. Therefore, the measurement of the vortex positions in different layers samples the same distribution and 
provides a signature for the expected transition into a squeezed state. In Fig. |5] we plot the relative squeezing as 
a function of rotation frequency at fixed temperature. As can be seen from this figure the size of the uncertainty 
ellipse can easily be several times the radial trap length and also becomes very strongly deformed as the rotation 
frequency is increased. Furthermore, the long axis of the uncertainty ellipse can be much larger than the vortex core 
size. The results in Fig. Elwere obtained with only one set of typical parameters. By a careful choice of, for example, 
lattice depth, number of particles, number of layers, or temperature a considerable degree of tunability is possible. 
For example, a shallower lattice implies a larger relative squeezing. 



IV. PHASE FLUCTUATIONS 



For our theory to be valid, constraints have to be set for the temperature. Since we assume a pure condensate, we 
are ignoring the influence of the noncondensate atoms on the vortex line. Therefore, the temperature should be much 
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FIG. 1: Solution of the gap equation as a function of temperature and rotation frequency, (a) Shows the order parameter as a 
function of temperature for two different rotation frequencies. The solid line represents il — Qc = 0.01 ■ cu,- and the dashed line 
represents fl — ^lc = 0.02 ■ cu,.- The result for the order parameter is scaled to the precession frequency of the vortex, (b) Shows 
the order parameter as a function of rotation frequency at two different temperatures. The solid line represents a temperature 
T = 10 nK and the dashed line represents the temperature T — 20 nK. The curves are calculated for the parameters A = 795 nm, 
LOr = 2iT ■ 100 Hz, A'^ = 1000, TVs ~ 51, and Vl = 15 -Br, where Er = /mX^ is the recoil energy of the atom after absorption 
of a photon from the laser beam. These parameters are representative of current experimental capabilities. We use this same 
set of parameters also in the other figures. 



lower than the critical temperature for Bose-Einstein condensation in the layers. This condition is relatively easy to 
satisfy. However, a more stringent condition is set by the phase fluctuations of the order parameter. In an infinite 
one-dimensional system phase fluctuations destroy the long-range order. In the finite system we are considering here, 
phase fluctuations can only be excited if the temperature is high enough psL l29l l30l | . In order to calculate the energy 
cost for a phase gradient of the order parameter we must determine the associated stiffness or superfluid density 
ps{T). Since the lattice breaks the Galilean invariancc this calculation is not entirely standard. 

The superfluid density is defined by the systems response to an imposed phase gradient '3l1. Therefore, we want 



to study the long-wavelength effects of the transformation v„ 
nearest-neighbor coupling and result in an interaction energy 



Vne 



i4>„/2 



The phase gradients enter only in the 



Jv 
2 



(n.m) 



0m) /2 - (0„ - (t)mY 



which can be considered as a disturbance in the long- wavelength limit. With this interaction energy we calculate in 
second-order perturbation theory the contribution to the free energy of the phase gradients. Up to second order the 
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FIG. 2: (Color online) Squeezing of the quantum mechanical position uncertainty ellipse. The figure shows the relative squeezing 
e = {y^ — i^)/ {x^ + jf) of the uncertainty ellipse of the vortex position as a function of rotation frequency at a low temperature 
of 5nK. Subplots show samples of the Gaussian probab ility dis tributions for the vortex position for two different rotation 
frequencies. In the subplots we use the trap length Ir = yhjmjjjr as a unit of length. In order to ease visual inspection we use 
different aspect ratios for the axes of the subplots. For the parameters used the coherence length at the center of the condensate 
is about 0.34 /r- 



effective imaginary time action for the phase fluctuations is 

= dT{Hj (r)) - - drdr'iHj (r) Hj (r')). 

Furthermore, it is sufficient to keep only the terms up to second order in the phase differences between layers. In this 
way we find the action for the phase gradients as 

rfr-^^ {cos(fc'A/2) [M^Nk' + \vk'\^ {Nk- + 1)] - JvPsin^ (fc'A/2) TVfc, {Nk^ + 1)} x ^fc2^fc</)_fc. 

■^0 k' k 

_ (11) 

This result should be compared with the standard (Galilean invariant) Landau |27j result 

Sl = dr j dz^p,{T)\V(^\\ (12) 
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where the superfluid density at temperature T is 



k 



\uk?Nk + bfep {Nk + 1) - —k^Nk {Nk + 1) 



(13) 



When we take the lattice spacing to zero while keeping Jy\^ constant, we can see that Eq. recovers the standard 
result when the mass of the pairing field is defined as ttia — ^tjik, where 

is the efi^ective mass of the kelvons. Furthermore, from this we see that the superfluid density in our case is defined 
by 

PsiT) = {cos(fcA/2) [\uk\^Nk + \vk\^ {Nk + 1)] - Jy/Ssin^ {kX/2) Nk {Nk + 1)} . (15) 

* k 

Equating the energy cost due to a An/NgX phase gradient with the thermal energy gives us an estimate for the 
temperature scale T,p = JvpsiT)'!r^X/2NskB of the phase fluctuations. In order to avoid phase fluctuations, we 
should be well below this temperature scale. 

At temperatures higher than the vortex can still be squeezed, but the phase fluctuations result in different main 
axis of the vortex position uncertainty ellipse in different layers. In Fig.|21we plot the critical temperature for Bosc- 
Einstein condensation, the critical temperature for the transition into a squeezed vortex, and the temperature scale 
of the phase fluctuations for typical parameter values. This figure indicates that the vortex spontaneously squeezes at 
a temperature that is easily accessible experimentally. Below this temperature there is a region in which the vortex is 
squeezed, but with fluctuating main axes. At even lower temperatures the phase fluctuations become negligible and 
long-range order over the size of the system is established. 



V. SUMMARY AND CONCLUSIONS 



We studied the equilibrium squeezing of the vortex line in an optical lattice and predicted that strong squeezing 
is indeed possible in the experimentally realistic parameter regime. Although the kelvon interaction induced by the 
vortex displacement is the most important one, other mechanisms for kelvon interactions also exist. In principle, 
the kelvons are also coupled to the collective modes of the Bose-Einstein condensate, for example, to the quadrupole 
mode 24, 32]. Such processes typically induce an effective attractive kelvon interaction and thus renormalize the 
value of the interaction strength, but do not make it positive. Therefore, we expect that the squeezing transition is 
robust with respect to coupling to the collective modes of the condensate. 



Acknowledgments 

This work is supported by the Stichting voor Fundamenteel Onderzoek der Materie (FOM) and by the Nederlandse 
Organisatie voor Wetenschaplijk Onderzoek (NWO). 



[1] R. J. Donnelly, Quantized vortices in Helium II (CUP, Cambridge, 1991). 

[2] G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. Mod. Phys. 66, 1125 (1994). 
[3] M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 83, 2498 
(1999). 

[4] K. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. Rev. Lett. 84, 806 (2000). 
[5] J. R. Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ketterle, Science 292, 476 (2001). 

[6] E. Hodby, G. Hechenblaikner, S. A. Hopkins, O. M. Marago, and C. J. Foot, Phys. Rev. Lett. 88, 010405 (2002). 
[7] D. A. Butts and D. S. Rokshar, Nature 397, 327 (1999). 
[8] J. E. WiUiams and M. J. Holland, Nature 401, 568 (1999). 

[9] R. J. Dodd, K. Burnett, M. Edwards, and C. W. Clark, Phys. Rev. A 56, 587 (1997). 
[10] F. ZambeUi and S. Stringari, Phys. Rev. Lett. 81, 1754 (1998). 



7 




FIG. 3: The phase diagram for the vortex squeezing transition. The upper most sohd hne is the critical temperature for the 
squeezing transition. The assumption of a pure condensate implies temperatures well below the critical temperature of the 
Bose-Einstein condensation indicated by the dashed line. Phase fluctuations can be ignored well below the lowest line. 



[11] A. L. Fetter and A. A. Svidzinsky, J. Phys.iCondens. Matter 13, R135 (2001), and references therein. 
[12] F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, Science 293, 
843 (2001). 

[13] S. Burger, F. S. Cataliotti, C. Fort, F. Minardi, M. Inguscio, M. L. Chiofalo, and M. P. Tosi, Phys. Rev. Lett. 86, 4447 
(2001). 

[14] C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, and M. A. Kasevich, Science 291, 2386 (2001). 
[15] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. Bloch, Nature 415, 39 (2002). 

[16] M. Grcincr, O. Mandel, T. W. Hansch, and I. Bloeh, Nature 419, 51 (2002). 

[17] D. Jaksch, C. Brudcr, ,]. Cirac, C. W. Gardiner, and P. ZoUer, Phys. Rev. Lett. 81, 3108 (1998). 

[18] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. A 63, 053601 (2001). 

[19] E. Cornell, private communication. 

[20] A. Fetter, Phys. Rev. 162, 143 (1967). 

[21] A. Fetter, Phys. Rev. 163, 390 (1967). 

[22] V. Bretin, P. Rosenbusch, F. Chevy, G. V. Shlyapnikov, and J. Dalibard, Phys. Rev. Lett. 90, 100403 (2003). 
[23] T. Mizushima, M. Ichioka, and K. Machida, Phys. Rev. Lett. 90, 180401 (2003). 
[24] J.-P. Martikainen and H. T. C. Stoof, Phys. Rev. Lett. 91, 240403 (2003). 

[25] Alternatively, the incomplete gamma function can be expressed in terms of the perhaps more familiar exponential integral 

as r[l -p,z] = z^-^Epiz). 
[26] H. T. C. Stoof, Phys. Rev. A 49, 3824 (1994). 

[27] J. R. Schrieffer, Theory of Superfluidity (Benjamin, New- York, 1964). 

[28] D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 85, 3745 (2000). 
[29] J. O. Andersen, U. A. Khawaja, and H. T. C. Stoof, Phys. Rev. Lett. 88, 070407 (2002). 



8 



[30] U. A. Khawaja, J. O. Andersen, N. P. Proukakis, and H. T. C. Stoof, Phys. Rev. A 66, 013615 (2002). 
[31] A. M. Rey, K. Burnett, R. Roth, M. Edwards, C. J. Williams, and C. W. Clark, J. Phys. B 36, 825 (2003). 
[32] J.-P. Martikainen and H. T. C. Stoof, cond-mat/0402052, [accepted for publication in Phys. Rev. A]. 



